A Lemaitre— Tolman— Bondi cosmological wormhole 
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We present a new analytical solution of the Einstein field equations describing a wormhole shell of 
zero thickness joining two Lemaltre-Tolman-Bondi universes, with no radial accretion. The material 
on the shell satisfies the energy conditions and, at late times, the shell becomes comoving with the 
dust-dominated cosmic substratum. 
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INTRODUCTION 

Static and asymptotically flat wormhole solutions of 
the Einstein equations have been known for a long time 
[l| . The study of wormholes developed with the seminal 
paper by Morris and Thorne [2|, after which many solu- 
tions were discovered (see [3j for an extensive discussion). 
The possibility that inflation in the early universe may 
enlarge a Planck size wormhole_to a macroscopic size ob- 
ject was contemplated in Ref. [J|. Dynamical wormholes 
were discovered and studied in Refs. [5|, |6[ and worm- 
holes in cosmological settings were contemplated in var- 
ious works ([7|, [8( and references therein), with particu- 
lar attention being paid to wormholes with cosmological 
constant A, which are asymptotically de Sitter or anti-de 
Sitter according to the sign of A [9( . 

After the 1998 discovery of the present acceleration of 
the universe [10( and the introduction of dark energy in 
cosmological theories to account for this cosmic accel- 
eration, there were claims that phantom energy, an ex- 
tremely exotic form of dark energy with P < —p (where 
P and p are the pressure and the energy density, respec- 
tively) could cause the universe to end with a Big Rip sin- 
gularity at a finite time in the future [11| . There was then 
a claim in the literature [12[ that, if a wormhole accretes 
phantom energy, it grows to enormous size faster than 
the background universe, swallowing the entire cosmos 
which would then tunnel through the wormhole throat 
and re-appear in a different portion of the multiverse be- 
fore reaching the Big Rip singularity. This claim was 
based on qualitative arguments and was later disproved 
by two classes of exact solutions of the Einstein equa- 
tions representing wormholes embedded in a cosmological 
background dominated by phantom energy [13J. These 
wormholes accrete phantom energy but, even if their ex- 
pansion rate differs from that of the cosmic substratum 
initially, they become comoving with it as the scale factor 
of the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
universe in which they are embedded grows. 

The first class of solutions consists of a zero-thickness 



shell which carries exotic energy and does not perturb 
the two copies of the FLRW universe which it joins. The 
second, more realistic, class is described by a general- 
ized McVittie metric [14| with an imperfect fluid and a 
radial energy flow, with the mass of the wormhole shell 
distorting the surrounding FLRW metric [13J. Another, 
less general, solution of the Einstein equations describing 
a cosmological wormhole comoving with the background 
was presented in Ref. [15I J . 

Cosmological wormholes are truly dynamical and inter- 
est in this kind of solution has developed in parallel with 
the increasing attention paid to cosmological black holes 
[16j . Additionally, gravitational lensing by wormholes 
was studied in 17| and numerical solutions interpreted 
as wormholes in accelerating FLRW universes were pre- 
sented in Refs. [18j. Recently, Maeda, Harada, and Carr 
have given precise definitions for general cosmological 
wormholes and have found two new exact solutions of 
this kind [19j. An important result of this work, which 
echoes a previous result of [7( , is that the null energy con- 
dition needs not be violated in this dynamical situation, 
although it must be violated for static wormholes to exist 
[19j |. It seems that the study of cosmological wormholes 
is developing into a promising new area of research. 

In this paper we propose a new analytical solution 
of the Einstein field equations describing a cosmologi- 
cal wormhole shell joining two Lemaitre-Tolman-Bondi 
(LTB) universes. We are led to this solution by the fol- 
lowing considerations: the second class of solutions in 
Ref. [13j is inspired by the McVittie metric, which de- 
scribes a central inhomogeneity in a FLRW background. 
However, the McVittie metric needed to be general- 
ized by removing the McVittie "no accretion" condition 
G*oi = (in spherical coordinates) which forbids radial 
energy flow. The goal of Ref. [13] was to describe the 
effect of the accretion of phantom energy onto the worm- 
hole. Here, we begin by noting that inhomogeneities em- 
bedded in a FLRW background are usually described by 
using an LTB metric [2fJ-|22[, not a McVittie one. The 
classical LTB metric describes a spherically symmetric 



inhomogeneity in a dust-dominated FLRW background. 
The Bondi condition Goi = parallels the McVittie no- 
accretion condition and forbids the (radial) flow of cosmic 
dust onto the inhomogeneity. It would be interesting to 
obtain a solution describing a wormhole shell joining two 
identical LTB universes and perturbing its surroundings 
in the way described by the LTB metric. This is what we 
do here. We obtain a wormhole shell composed of exotic 
matter which expands more slowly than the cosmic sub- 
stratum (which becomes a spatially flat FLRW universe 
at late times), but eventually becomes comoving with it. 
The next section details how to construct the wormhole 
shell and satisfy the Israel-Darmois-Lichnerowicz junc- 
tion conditions |23| on this shell. The Einstein equations 
on the shell provide expressions for the energy density 
and pressure of the material on the shell. Sec. Ill uses 
the covariant conservation equation to relate the rate of 
change of the mass of shell material, the shell area, and 
the flux of cosmic fluid onto the shell due to the relative 
velocity between the shell and the cosmic substratum. 
The metric signature is — \- -\ — \-, we use units in which 
the speed of light and Newton's constant are unity, and 
we follow the notations of Ref. [25|. Greek indices run 
from to 3 and Latin indices assume the values 0, 1, and 
2 corresponding to the coordinates (t, 9, <p) of the spher- 
ical hypersurface E defined below. 



THE LTB WORMHOLE SOLUTION 

The spherically symmetric LTB line element for the 
critically open universe in polar coordinates (t, r, $, ip) is 

a , r-D'n „m 2 j„2 , D 2 



ds 2 = -dt 2 + [R'(t, r)} dr 2 + R? (t, r) dtt 2 , (1) 



where 



R(t,r)=(rV 2 + ^^J7)?j 



is an areal radius, r is a comoving radius, 



m e {r) = An I dxx po(x) 



2/3 



(2) 



Pair) is the energy density on an initial hypersurface, 
a prime denotes differentiation with respect to r, and 
dVt 2 = dtf 2 + sin 2 •& dip 2 . The line element flTJ) describes 
a spherical inhomogeneity in a dust-dominated universe 
( 20h22|: for a recent review see |24j). 

Consider now a wormhole shell E at r = rj;(t) which 
joins two identical copies of an LTB spacetime (this shell 
describes a wormhole created with the universe and not 
formed as the result of a dynamical process after the 
Big Bang). The wormhole shell is dynamical and moves, 
possibly also relative to the cosmic substratum, and its 



motion is described by the form of the function r-s(t). It 
is convenient to write the equation of this shell as [27j 

f(t,r) = R-R I1 (t,r I1 (t))=0. (4) 

To find the unit normal to £ we first compute 



N„ 



and 



V M / = V„ (R - Ife) = (R t - R s , R' , 0, o) , (5) 



N» = [R s -Rt, — ,0,0 



(6) 



aNjj,. Here R t 



and then normalize according to n^ 

dR/dt and an overdot denotes a total derivative with 

respect to t, i.e., R = dR/dt. 

The normalization n^n^ = 1 yields 



1 



1 - [Rt - 



Re)' 



(7) 



It is convenient to introduce the radial velocity of the 
wormhole shell relative to the cosmic substratum 



v= Ry 



R 



(S 



where R t , s = R t (t,r E (t)). Then 



a 



VT 



liv) 



(8) 



(9) 



is an (instantaneous |28| ) Lorentz factor for the relative 
motion shell-background. The unit normal is then 

n^ = (-7«,7JZ',0,0), (10) 

n" = (7^,0,0) . (11) 

The restriction of the metric to E is given by 

dsL = -dt 2 + R'^drh + R^dQ 2 (12) 



(3) or, using the fact that R-z = i? t i s + i? E ^s on the shell, 



ds h = -(l-v 2 )dt 



2 ^ u2 ' Rldn 2 . 



(13) 



which expresses the fact that the proper time r of the 
shell is given by 



dr = v 1 — v 2 dt, 



(14) 



i.e., it is Lorentz-dilated with respect to the comoving 
time t of the background. 
Using the triad 

{e«. e W' e M } = {VT^ 2 S at ,6 ai, ,6^} , (15) 



the extrinsic curvature of the shell is given by (a, b, c = 

t,d , or if) 

K aP = e<? ef V a n b = eg* eg (cU 6 - T c ab n c ) , (16) 

where T c ab are the Christoffcl symbols of the 3- 
dimensional metric g ab \s- Eq. (fl6|) yields 



(17) 



#tt = £ (ftn t - Tint) = -a. - 2 7Wf « 2 , 



if MM = ^,9,9 sin 2 $ . 



The mixed components of the extrinsic curvature are 

K\ = 7 « t (2 7 V + l) , (18) 



K $ = J V — — = A ^ 



V ) 



(19) 



while the trace is 



' : r."' VP -9-v 3 , ' * |E 



i?S 



+ Vtv + 7f t 



(20) 

Since there are two identical LTB universes joining at the 
shell with unit normal n^ pointing outward, the jumps 
of these quantities on S are 

[K a b ] = 2K a b , [K]=2K. (21) 

The Einstein equations at the shell £ are [26[ 

[K\-5\K] = -^S a b , (22) 

where S ab is the energy-momentum tensor of the mate- 
rial on the shell. We assume that this matter is a perfect 
fluid, described by 

Sab = {v + P*)uf)uf ) +Pv9ab\ll, (23) 

where a and Pe are the 2-dimensional surface density 
and pressure, respectively, and u^ is the 4-velocity of 
the shell given by 

a _ dx< (V) _ ® X (E) dx^ dx^ ^(E) 



(s) dr dx* 1 dr 



= 7 



dt dx» 



(24) 



The coordinates on £ are x 1 ^ = (i, rs(t), •&, ip), which 
yield 

',7^,0,0), u^ = (-7,7^^,0,0). 

(25) 
As such, it is easy to see that 



*(=) = T 



U (£) U f ] =1 > W (£) »/* = ° ' 



(26) 



The (t, i) component of the Einstein equations (1221) at 
the shell is 



3 v -Rtis 

while the ($, 1?) or the (y>, <^) component yields 

D T /" , o 2 2 ,2 R t\s 

Pt, = — I ft + 27 w w t + 7 «-^— 

a v t 1 + v 2 
~2 +7 4^1-*u 2 ' 



(27) 



Using eq. $2$, one obtains 



Ri 



V m e ( r ) 



# r 3 /2 -1- 3Vm7W t/2 ' 



(28) 



(29) 



this quantity is asymptotic to Jr, the Hubble parameter 
of the dust-dominated cosmological background, as t — > 
+00. It is also 



2^"Rs~ 






s/m e (r E ) 



2 ^ 4 /2 + 3v/m e (re) t/2 



(30) 



and a > is equivalent to v < 0. A wormhole shell 
with positive surface energy density must necessarily ex- 
pand slower than the cosmic substratum, a fact that is 
interpreted as the influence of the inhomogeneity slow- 
ing down the expansion locally. Since the expansion 
rate of the background 3r tends to zero at late times, 
the shell expansion rate must also tend to zero and the 
shell becomes comoving. In fact, since v < 0, it is 
Rz = R t \T, + R's^s < Rt\s and, since i? E > 0, it is 
rs < 0. This inequality is consistent, of course, with the 
relation 



r-n 



i?v 



which is easy to derive. 

Now, rs > decreases monotonically implying that 
29], as t — > +00, either r^ tends to a horizontal asymp- 
tote Too > 0, or lim t ^ +oc rs(t) — + . If r% —¥ + , then 



i?E 



(r¥ + \V^^)t) ^0 (31) 



because m e (0) = and the wormhole disappears asymp- 
totically, which doesn't make sense physically, and this 
possibility is discarded. Hence, 



Re 



(r% 2 + \V^t^o)t) (32) 



and the wormhole shell becomes comoving at late times. 
We conclude this section with a comment on the en- 
ergy conditions. The strong energy condition for a 2- 
dimensional perfect fluid is a + P^ > and a + 2Pe > 0: 



for our wormhole, assuming a non-negative surface den- 
sity a and hence v < 0, it is 



a "/v t 1 + v 



<? + Pe = -t + ^- 



2 4tt 1-v 2 



>0 



and 



7Vt 1 + V 2 

47T 1 — V 1 



(33) 



(34) 



because v < and w — > 0~ as £ — > oo, hence «t > 0. 
The weak energy condition on the shell corresponds to 
a > and a + Py, > 0, while the null energy condition 
is equivalent to a + P^ > 0. Therefore, the material 
on the shell satisfies the weak, strong, and null energy 
conditions. 



THE COVARIANT CONSERVATION EQUATION 



This formula is interpreted physically as follows: the 
quantity p v is the flux density of cosmic fluid onto the 
shell caused by the relative motion of the shell with re- 
spect to the background. The quantity pvA^ is the flux 
of this material; the factor 2 appears because there are 
two LTB spacctimes joining at the shell. The Lorentz 
factor 7 is due to the Lorentz contraction caused by the 
radial motion of the shell. 

Eq. (|4"Tj) expresses the first law of thermodynamics re- 
lating changes over a time interval dt 



dM + dW = dg s 



(42) 



where dM — Mdt is the variation of internal energy dur- 
ing dt, dW = Ps A^dt is a work term due to the variation 
of the area of the shell, and dQs is the energy input due 
to the influx of cosmic fluid onto the shell. 



We can now solve the covariant conservation equation 
projected along the 4-velocity of the shell u?^ [26| , 



It is convenient to note that 



a rji 13 
u (T,) T c n P 



?r\VhS„ = M?v>V 



'(E) v & J a 



(E) v b 



C=)«L 



(<7 + P E )u^ E) 



(35) 



ufaVaPx 



-V6(<7U( E) ) -P S V 6 ^ E) 



(36) 



using the normalization «?„>«„ = — 1 and its conse- 



,( s 



'(£)• 



quence u^^V^Wq = 0. We now compute 



"(E) 



Vb ( au 



4 (S) 



5|E v 



ff|E|^M (E) 



(37) 



where ffm = 7 2 R% sin 2 fl is the determinant of the 3- 
dimensional metric g a b\Y,i obtaining 

V&(o-W(E)) - 



Rl ^l°)=l A 



M 



(38) 



Here A^ = 47ri?| is the area of the shell and M = aA% is 
the mass of the material located on the shell. Similarly, 
one obtains 



V7 b 

VbW (s) = 7 



As 

As 



(39) 



and 



rj a rp j3 



«(E)T„% =2nf E) T/n^ = -2 7 >. (40) 



Putting everything together, we obtain the covariant con- 
servation equation in the form 



M + Pe^s = -2 7 pt;A s . 



(41) 



DISCUSSION AND CONCLUSIONS 

We have obtained, and interpreted physically, an ex- 
act solution of the Einstein field equations representing 
a wormhole shell joining two identical LTB spacetimes 
which are dust-dominated. This solution is similar to 
the wormhole solution of Ref. [13j obtained by general- 
izing the McVittie metric, but there are important dif- 
ferences. First, we adopted the no-accretion condition 
of Bondi 22] which forbids radial flow of energy into the 
wormhole while Ref. [13|, being interested in the effect 
of accretion, allows for radial flow with the consequence 
that an imperfect fluid is needed in order to obtain so- 
lutions in 13]. Here, instead, we can c onsider a perfect 



fluid, the dust of classical LTB solutions [2CH22| . While in 
[la ] the conservation equation analogous to our eq. (f4"Tj) 
has a right hand side consisting of two terms, one due 
to the relative motion between shell and cosmic substra- 
tum, and another due to accretion, only the first term 
appears in our case in which there is no radial flux. 

An important result of [19( is that, contrary to static 
wormholes, the null energy condition needs not be vi- 
olated for their cosmological and dynamical wormholes 
to stay open; here we propose a different cosmological 
wormhole solution made with material which satisfies the 
weak, null, and strong energy conditions on the shell. In 
other words, the "stuff" necessary to keep this worm- 
hole throat open does not need to be very exotic. This 
feature motivates further studies of dynamical wormhole 
solutions of the Einstein equations. 
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